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This paper extends work done to date on quantum com- 
putation by association of potentials with different types of 
steps. Quantum Turing machine Hamiltonians, generalized 
to include potentials, correspond to sums over tight bind- 
ing Hamiltonians each with a different potential distribution. 
Which distribution applies is determined by the intitial state. 
An example, which enumerates the integers in succession as 
binary strings, is analyzed. It is seen that for some initial 
states, the potential distributions have quasicrystalline prop- 
erties and are similar to a substitution sequence. 

89.70.+c,71.23Ft 



I. INTRODUCTION 

Quantum computation is a field of much interest. Since 
the early work , much impetus was provided by 
Shor's discovery |H] that for some problems quantum 
computers are more efficient than classical computers. 
Recent work includes discussion of quantum error correc- 
tion Q, use of quantum gates for computation ||], and 
Hamiltonian descriptions of quantum Turing machines 




In all work to date models of quantum computation 
are considered in which successive steps of the compu- 
tation proceed smoothly with no potentials present, if 
the system is isolated from the environment (see for ex- 
ample 1^). Landauer |^ has noted that environmental 
influences and defects in physical models introduce ran- 
dom potentials which degrade performance by causing 
reflections at various steps and decay of the transmitted 
component. 

Here potentials are introduced into quantum compu- 
tation, not in a random fashion, but by association with 
different types of steps. An example is the association of 
a potential with all Turing machine "read 1" steps. All 
other steps are potential free. 

The method used is to generalize the description of 
quantum Turing machines to include potentials at 
one or more types of computation steps. Isolation from 
the environment is assumed. In this case the Feynman 
Hamiltonian for these generalized quantum Turing 
machines (GQTM)s, can be also be described as a sum of 
tight binding Hamiltonians each of which has a different 
potential distribution along a one dimensional path of 



computation states. Which of these Hamiltonians applies 
to the GQTM is determined by the initial state. 

For some GQTMs the potential distributions are sim- 
ilar to those which are of much recent interest in con- 
densed matter physics. This includes Hamiltonians with 
quasiperiodic potentials jl^ or with potential distribu- 
tions corresponding to substitution sequences (see |Tll ] 
for a definition). Examples include the Period Doubling, 
||l|,|3|, Thue- Morse (l|Jl|l, Fibonacci @, and Rudin- 
Shapiro IT^ ] substitution sequences. This similarity pro- 
vides a link between the behavior of a GQTM as it evolves 
during a computation and the 1-D motion of a particle 
such as an electron in systems with potentials which are 
quasiperiodic or correspond to substitution sequences. 

This similarity will be shown here by first giving a 
brief description of the physical model and the sum de- 
composition of the GQTM Hamiltonian. An example, 
the counting GQTM which enumerates the nonnegative 
integers in succession as binary strings will be considered. 
This is of interest because some of the potential distri- 
butions are similar to but not identical with those for 
both quasicrystals and substitution sequences. The dif- 
ferences suggest that the distributions are a new type of 
1-D structure that does not seem to have been discussed 
in the literature. 

II. THE PHYSICAL MODEL 

For quantum Turing machines, the physical model cor- 
responds to one tape machines. It includes a finite state 
head that moves along an infinite lattice of finite state 
systems called qubits. For many computations it is suf- 
ficient to consider binary strings corresponding to qubit 
states |0), |1). The actual physical form of each qubit 
(e.g. as a spin 1/2 system, etc.) is not important here. 

The computation basis B which spans a separable 
Hilbert space 7i is the set {\l,j,S)} of states. Here j 
and I denote the lattice location and internal state of the 
head. 15*) = (^J^-ool^U)) denotes the lattice qubit state 
where S is any function from the integers into the set of 
possible qubit states such that S{j) ^ for at most a 
finite number of j values. This condition, the tail state 
condition, is one of many that can be imposed to keep 
B denumerable. Both B and Ti depend on the condition 
imposed. 



1 



III. STEP OPERATORS AND HAMILTONIANS 

Each GQTM is described by a step operator T such 
that iteration of T or its adjoint corresponds to successive 
steps of the GQTM in the forward or backward time 
direction. T is defined as a finite sum of elementary step 
operators, 7" = X); s = J2is H,sWi^s, over head states 
I and qubit states s where ji^s is a real constant between 
and 1.. Wi^s has the form, 

oo 

W;,s = ^ WusQlVLsPs.jUi^sPj (1) 
j = -oo 

where Q;, Pj, Pgj are respective projection operators for 
the head in internal state \l), at site j, and the site j 
qubit in state \s). 

The unitary operators wi^sjVi^sjUi.s describe head and 
qubit state changes, and head motion. The depen- 
dence on I, s is indicated by the subscripts. The op- 
erators wi^s and ui^s satisiy the commutation relations, 
wi^sQi = Qi'WLs for some I'eL and ui^gP-j = Pj±iui^s- 
Here L is the finite set of head states and Pj+l^ Pj^i cor- 
respond to head motion one site to the right or one site 
to the left. For our purposes it is sufficient to require 
that vi^s satisfy vi^gPsj — Ps'jvi^s where s'eS, the set of 
qubit states. If the qubits are binary this is equivalent to 
requiring that u;.s = or vi,s = 1. 

These step operators generalize earlier definitions 
in that values of 7;^^ different from or 1 are al- 
lowed. Potentials are associated with steps for which 
terms Ti^^ with 7;^^ < 1 are active. 

The main condition imposed on T is that it be distinct 
path generating on the computation basis B. This condi- 
tion can be relaxed to apply to any basis, not just B (see 

for details) but this generality is not needed here. The 
condition means that iteration of T or on the states 
in B generates finite or infinite paths of states which do 
not join, branch, or intersect. Mathematically this is ex- 
pressed by requiring that T be a direct sum of weighted 
shifts That is T = UD = UzD, where for each 

1 Ui is a bilateral shift, a unilateral shift, the adjoint of 
a unilateral shift, a finite truncated shift of length TV, 
or a cyclic shift of length M with N, M arbitrary. D is 
defined by D\l,j,S) = 7;,S0)l^j, S) for any state \l,j, S) 
in B not annihilated by U or W . Both U and D satisfy 
U = J2i UiPi,D = X^ DiPi where the Pi are orthogonal 
multidimensional projection operators on H. 

The i sum is a sum over paths where each path i in 
B is defined by iteration of Ui or . That is if the 
state \l,j,S) is not annihilated by Ui or Uj then path i 
is defined by ■■■ ujH, j, S) ,\l, j, S),U\l, j, S) ,■■ 

For each step operator T define a Hamiltonian by , 

H = K{2-T -T'') (2) 

where K is an arbitrary constant. If T = UD, then 
H can also be written as if = J^i HiPi where Hi — 



K{2 — Ti — T/) describes motion along the path states in 
Bi. Eq.|| is based on the association of an infinitesimal 
time interval with T so that it can be directly used to 
construct a Hamiltonian. 

From this it follows that for each path i, Hi is a tight 
binding Hamiltonian on PiTC with nearest neighbor off 
diagonal potentials. To see this note that 

H,^K{2-U^~Ul) + V, (3) 

where the off-diagonal potential Vi along path i is given 
by V, = K[U^{1 - A) + (1 - D,)U}]. The constant term 
2K is present so that for each path, the path kinetic 
energy component of Hi, K{2 — Ui — Uj) is equivalent to 
the symmetrized discrete version of the second derivative, 
~K(f/{dxf. 

Under the action of H the initial state ^'(0) deter- 
mines which tight binding Hamiltonian applies. If ^'(0) 
is a linear sum of computation states in PiU., then Hi 
gives the evolution of ^1/(0). If \E'(0) has components in 
more than one subspace PiTi, then each component has 
a corresponding Hi which describes its evolution. 

For many i the potential distributions in the corre- 
sponding Hi are the same. This is the case for any 
transformation group whose actions do not change T. 
For example, for GQTMs, each shift of an initial state 
along the lattice corresponds to a different path i with 
the same Hi . From this one sees that equivalence classes 
Vi = {Pj : Hj = Hi} of projectors can be defined which 
identify paths with the same Hi. 

It is important to distinguish between path motion and 
head motion on the qubit lattice. In general, iterations 
of T describe head motion backward and forward over 
lattice sites. However, (cyclic paths excepted) for each 
passage over a site, the overall system state is orthogonal 
to the state for any other passage. As a result the back 
and forth lattice motion can be unfolded to motion along 
a path where the path states are pairwise orthogonal and 
are ordered by iteration of T or . The tight binding 
Hamiltonians, Eq. |^, refer to motion along paths and not 
the qubit lattice. 

IV. THE COUNTING GQTM 

Simple examples can be constructed by letting 7;^^ as- 
sume just two values, 1 and 7 where 7 is a real number 
between and 1. In particular, for all leL, let 7;_i = 7 
and 7;^s = 1 for all s ^ 1. 

This simplification corresponds to GQTMs in which a 
potential is associated with all "read 1" elementary steps 
(those steps with P^j = P\j in Eq. Q). All other qubit 
state read actions are potential free. In this case the I, s 
sum for T can be split into two parts, T = Wi s + 

The counting GQTM can be used to illustrate the fore- 
going. This GQTM, starting with a blank qubit lattice 
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(all Os) with two markers spaced n+1 sites apart, gener- 
ates all binary numbers in succession from up to 2" — 1 
by continually adding 1. 

The step operator for this GQTM, which is used in Eq. 
H, is a sum of 7 terms: 

oc 
j=-oo 

+ wQiP2jU^Pj + ^Q2VxjPljU^ Pj + w'^Q2VxjPojUPj 
+ wQ2P2jUPj) (4) 

The projection operators are as defined for Eq. [|; w 
is a shift mod 3 on the three head states {wQm — 
Qrn+iw mod 3) and u shifts the head along the lattice 
by one site {uPj = Pj^iu). The need for markers is 
accounted for here by choosing the qubits in the lat- 
tice to be ternary with states |0),|1),|2). |2) is used 
as a marker and |0}, |1) are used for binary strings. 
The qubit transformation operator for the site j qubit 
Vxj — o-xj{Poj + Pij) + P2j exchanges the states |0), |1) 
and leaves the state |2) alone. 

The terms in T are chosen so that T is distinct path 
generating on B. The theorems given in 0|, that verify 
this conclusion were proved for D ^ 1. The proofs should 
also hold for _D 7^ 1. 

An initial state for this GQTM is shown in Figure 1 
with the head in state |0} in a wave packet localized to 
the left of the origin. All qubits are in state |0) except 
those at sites 0,7i -f 1 which are in state |2}. Succes- 
sive iterations of T move the head to the righthand units 
marker followed by successive enumeration of the first 2" 
numbers as binary strings. When the space between the 
markers is filled with Is, the initial state is restored by 
conversion of the Is to Os. The head in state |1) moves 
to the right with no more qubit or head state changes. 

For this GQTM, a potential of height V = 2K{1 - 7) 
is present at path states at which term 5 is active. This 
is the only term with a "read 1" operator. The path 
potential distribution is determined by a close examina- 
tion of successive iterations of T. For example, to add 

1 to • • -01001112 to obtain • • -01010002 with the head 
beginning and ending in state |1) at the units marker 

2 requires one "read 2" step (term 4) followed by three 
"read 1" plus 1 — » steps (term 5) followed by one "read 
0" plus — > 1 step (term 6) followed by three "read 0" 
steps (term 3). In this path segment there is a potential 
of height V and width 3 with potential free regions of 
width 1 to the left and width 4 to the right. Extension 
of the segment shows a few more potential free sites in 
both directions before encountering the next potential. 

The path potential distribution is expressed in general 
as follows: Let i? be a function from the nonnegative 
integers to the nonnegative integers such that R{j) gives 
the number of Is occurring before the first in the binary 
string for j. R{j) = if j even and R{1) — 1, i?(3) — 



2, i?(5) = 1, R{7) = 3, - - - . R can be expressed recursively 

by 



R71 — 'S'n — 1?^; Sn — RnSn 



(5) 



for n — 1,2,-- - with 6*0 = where i?„ is the initial 
segment of R of length 2" . 

It turns out that i? is a substitution sequence pT| with 
the substitution rule n — > On + 1 for n = 0, 1, 2, - - -. (That 
is, the infinite sequence 0, 1, 0, 2, 0, 1, 0, 3, • - - is invariant 
under the replacement of n with On -f 1.) This gener- 
alizes the literature definition of substitution sequences 
|l3| , ^ 19 1 in that the alphabet is infinite. 

A study of the iteration of T with the initial head state 
1,71 + 1) and 15*) as shown in Fig. 1 shows that the 
potential distribution is given by a binary sequence TZ 
obtained from R by replacing R{j) by a string of numbers 
according to the following prescription: 00 = and 
n OnO"^"'^ = n. Here 0"+^ denotes a string of n -I- 1 
Os and n denotes a string of n Is. Justification for this 
can be seen by extension of the example of adding 1 to 
•--01001112. 

TZ is generated recursively by Eqs. [s] if n is replaced 
by n and Sq — 0. It follows that TZ is also a substitution 
sequence for the words 0, 1, • - - regarded as an alphabet. 
The recursive generation is important for future work on 
the spectral and transmission properties of tight binding 
Hamiltonians with this potential distribution. 

The above is illustrated in Figure 2 for n = 6. The fig- 
ure shows the distribution of read and read 2 steps (as 
Os) and read 1 steps (as Is) as a function of the number of 
T iterations starting from the initial head state |1, rt + 1) 
and the two markers separated by 6 sites (Figure 1). Fig. 
2 also gives the potential distribution for enumerating 
the first 64 numbers as binary strings. Changing the ini- 
tial state to any component of the wave packet shown in 
Figure 1 shifts the distribution along the path with no 
changes in the width or spacing between the potentials. 

Different initial states correspond to different poten- 
tial distributions. Changing the separation n + 1 be- 
tween the two markers changes the number of poten- 
tials present. Initial states with m markers separated 
by rii + 1, + 1, - - - J n-m + 1 qubit sites give a poten- 
tial distribution obtained by concatenating (with some 
intervening Os) initial segments of TZ corresponding to 
enumeration of the first 2"^ , 2"^ , . . . , 2"™ numbers. If all 
the rii — p are the same, the potential distribution is 
periodic with m periods each containing the first 2^"^ 
potentials. 

If just one marker is present and the head is in state 
|1), then all binary numbers are generated in succession 
without halting. For this input state the potential dis- 
tribution in Hi, Eq. ||, is given by TZ. The system is 
quasicrystalline in that no two path sites have the same 
global potential environment, yet each local environment 
of arbitrary length in TZ is repeated infinitely often pOl . 
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DifTcrcnccs between this distribution and those described 
in the hterature include the facts that the potential distri- 
bution is one-way infinite, the alphabet is infinite, and IZ 
is a substitution sequence with respect to alphabet words 
and not the original alphabet. Because of this the 1-D 
structure appears to be of a new type not yet analyzed. 
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FIGURE CAPTIONS 

Figure 1. Initial and Final States for Counting GQTM 
for the First 2" Binary Numbers. All lattice qubits are 
in state |0) except those at sites and n + l which are in 
state 1 2). The initial and final head states are shown as 
wave packets with internal head states ]0) and |1) to the 
left and right respectively. 

Figure 2. Path Potential Distribution for Enumeration 

of the first 64 Numbers (n = 6) as Binary Strings. The 
figure shows the distribution of Os, corresponding to read 
or read 2 steps, and Is, corresponding to read 1 steps, 
as a function of the number of T iterations starting from 
the initial head state |1, n -|- 1). 
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